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ABSTRACT 

We investigate the dynamics of self-gravity wakes in dense planetary rings. In par¬ 
ticular, we examine how the pitch angle of self-gravity wakes depend on ring parameters 
using A-body simulations. We calculate the pitch angles using the two-dimensional au¬ 
tocorrelation function of the ring surface density. We obtain the pitch angles for the 
inner and outer parts of the autocorrelation function separately. We confirm that the 
pitch angles are 15 to 30 degrees for reasonable ring parameters, which are consistent 
with previous studies. We find that the inner pitch angle increases with the Saturnicen- 
tric distance, while it barely depends on the optical depth and the restitution coefficient 
of ring particles. The increase of the inner pitch angle with the Saturnicentric distance 
is consistent with the observations of the A ring. The outer pitch angle does not have 
the clear dependence on any ring parameters and is about 10-15 degrees. This value is 
consistent with the pitch angle of spiral arms in collisionless systems. 

Subject headings: planets: rings, methods: A-body simulations 


1. Introduction 


The main rings of Saturn are composed of numerous icy particles and exhibit many different 
patterns, including non-axisymmetric structures on sub-kilometer scales. Gravitational forces be¬ 
tween particles form gravitationally bound clumps, while differential rotation tears them apart. Due 
to these competing processes, spatial structures called gravitational or self-gravity wakes appear 


(Colombo et al. 1976; Salo 1992a). Salo (1992a) performed the first investigations of self-gravity 


wakes using A-body simulations. The typical radial scale of the wakes is comparable to the critical 


wavelength of gravitational instability (Toomre 1964), which is around several tens to hundred 


meters. Their azimuthal scale is much longer than their radial scale. 
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Even though individual wakes are tiny structures, together they can have detectable influence 
on the large-scale photometric properties of the rings. The A and B rings exhibit a remarkable 


asymmetric brightness variation that was first discovered by Camichel (1958), confirmed later by 


the ground-based observations ( Ferrin||1975 Reitsema et al.||1976 ), and by the Voyager image data 
( ]Franklin et al. 1987; Dones et al. 1993). Salo & Karjalainen (2003) and Salo et al. (2004) carried 
out the detailed photometric modeling of self-gravity wakes with A^-body simulation. They found 
that the self-gravity wakes in A^-body simulation can account for the brightness dependence on 
the longitude and elevation in the Voyager and ground-based observations. This conclusion was 
confirmed by the Hubble Space Telescope and Cassini observations (French et al.|2007 Porco et al. 


2008). 


Though there have been many numerical A-body simulations of self-gravity wakes ( 

Salo|1992b 

Richardson 

1994 

Daisaka &: Ida 

1999 

Daisaka et al. 

2001 

Robbins et al. 

2010 

Yasui et al. 

2012), 


the theoretical framework for dense planetary rings with self-gravity wakes has not been explored 
in detail yet. It is often argued that wakes in dense rings are dynamically similar to the spiral 
structure in differentially rotating galactic disks (Colombo et al.|[l976 Franklin et al.|[T987). Julian 


&: Toomre (1966) studied theoretically the spiral structure in disk galaxies. They found that when 


disk surface density is perturbed, non-axisymmetric density enhancements appear. This mechanism 
is called swing amplihcation ( [Toomre 1981). Based on the swing amplification, we can estimate 
the pitch angle of a spiral arm ( Julian &: Toomi^|1966 Toomre||1981 ). The pitch angle is the angle 
between the tangents to a spiral arm and a perfect circle, which measures how tightly the spiral 
arm is wound. The swing amplification model predicts that the pitch angle for Keplerian rotation 
is about 10-15°, which is roughly consistent with that of the wakes in rings. Thus, it has been 
considered that the swing amplihcation plays an important role in the self-gravity wake formation 
in rings. 

According to the swing amplihcation model, the pitch angle of the self-gravity wakes depends 


on the shear rate (Julian &: Toomre 1966; Michikoshi Sz Kokubo 2014). However, the dependence 


of the pitch angle on other ring parameters is still not well understood. Observationally [Hedmaii 


et al. (2007) investigated the dependence of the pitch angle on the Saturnicentric distance. They 


found that the pitch angle varies systematically with the Saturnicentric distance. In the A ring. 


the pitch angle is 18°-24°, which is consistent with the previous numerical simulations (Salo et al 


2004). In the middle A ring, the pitch angle is almost constant 18°, but in the outer A ring, it 


increases with the Saturnicentric distances to 24°. In addition, Hedman et al. (2007) pointed out 
that the pitch angle has the local maximum at the location of the strong density waves. From 


the Cassini UVIS occultation studies, Colwell et al. (2007) obtained the consistent results that the 


pitch angle in the A ring is 10-50° and it increases with the Saturnicentric distance. They also 
found that the pitch angle in the B ring decreases with the Saturnicentric distance. From the HST 


observations, French et al. (2007) showed the increase of the pitch angle with the Saturnicentric 
distance in the A ring, which is consistent with the results in Colwell et al. (2007). The pitch 


angle dependence on the Saturnicentric distance has a complex structure. It is not a monotonic 




























































































function of the Saturnicentric distance. However, roughly speaking, the pitch angle decreases with 
the Saturnicentric distance in the B ring, while it increases in the A ring. 


Salo (1995) performed Al-body simulations and found that the wake pitch angle is about 20^ 


from the autocorrelation function. 


Salo et al. (2004) found that the pitch angle varies with the 

and 


autocorrelation distance. The pitch angle for small autocorrelation distance is about 25-30' 
decreases with the distance towards 15-20°. Using the photometric model they found that the 
longitude of the minimum brightness with respect to ansa increases from 17° to 23° with the eleva¬ 
tion angle of the observer with respect to the ring plane. They argued that the observed minimum 
longitude for high elevation angles corresponds to the pitch angle for small distance, and that for 
low elevation angles corresponds to the pitch angle for large distance since the brightness for high 
elevation angles is dominated by the inner dense part of the wakes. However, the systematic study 
of the dependence of the pitch angle on ring parameters has not been conducted yet. Therefore, 
we perform such a survey with high resolution A^-body resolutions. We calculate the pitch angle 
from the autocorrelation function of the ring surface density and investigate its dependence on the 
restitution coefficient of ring particles, the dynamical optical depth, and the Hill radius of a particle 
pair relative to the sum of their physical radii. 


The outline of this paper is as follows. We describe the simulation model in section 2 and 
explain the calculation method of the pitch angle in section 3. In section 4, we present the results 
of the dependence of the pitch angle on ring parameters. Sections 5 and 6 are devoted for discussion 
and a summary, respectively. 


Simulation Method 


2.1. Equation of Motion 


We adopt the shearing box method developed for local ring dynamics (e.g., Wisdom & Tremaine 
1988; Michikoshi Sz Kokubo|[2C)ll ). We introduce the rotating local Cartesian coordinates (x, y, z). 
The x-axis points radially outward, the y-axis points in the direction of orbital motion, and the 
z-axis points to the normal to the ring plane. We assume <C Ms, |xi|, |?/i|, |zi| <C o, and 
|dxj/dt|, |dyi/dt|, \dzi/dt\ <C all, where Ms is the Saturn mass, rrii is the mass of ith ring particle, 
a is the Saturnicentric distance, and H is the orbital frequency. Then the equation of motion for 
the ith particle is (e.g., [Petit &: Henon 1986, Nakazawa &: Ida||1988) 

























where is the distance between ith and jth particles. The first terms in Equations Q and Q 
are the Coriolis force. The second term in Equation 0 is the tidal force, the first term in Equation 
© is the vertical force due to Saturn’s gravity, and the last terms in Equations 0-@ are the 
mutual gravity between ring particles. 

The simulation box has a square shape with length L on the x-y plane. We use the sliding 
cell method where when a particle crosses a box boundary outward, the corresponding particle 


immediately comes into the box from the opposite boundary (Wisdom &: Tremaine 1988; Toomre 


&L Kalnajs [1991 ). The length L should be sufficiently larger than the characteristic scale length, 
the critical wavelength of the gravitational instability 


Acr — 


112 


( 4 ) 


where Eg is the mean surface density (Toomre 1964). We adopt L = 10A, 


Eor the collision model we adopt the constant restitution coefficient e that ranges from 0.1 to 
0.6. For e > 0.7, the velocity dispersion of ring particles increases monotonically and self-gravity 
wakes do not appear (Goldreich &: Tremaine||1978 ). The effect of the velocity-dependent restitution 


coefficient is discussed in Section 4.1 For the velocity-dependent restitution coefficient model, we 
use the standard model of Bridges et al.| (1984): 


e(un) = min(l, {vjvc) 


( 5 ) 


where is the normal component of the relative velocity in collisions and Vc = 0.0077 cms“^. We 
adopt the hard-sphere model as the collision model (e.g.. Wisdom & Tremaine 1988 Salo 1991} 
Richardson 1994). In the hard-sphere model, when collision between particles is detected, the 


particle velocities change instantly. In this study since we neglect the friction between particles, 
the relative tangential velocity is conserved in collision and the effects of spin are not considered. 

The equation of motion is integrated with a leapfrog integrator with timestep At = (27r/II)/200. 
We calculate the self-gravity by direct summation using the special-purpose computer, GRAPE-DR 
( Makino et al.||2007 ). We introduce the cutoff of self-gravity by the subregion method ( [Daisaka &: 
Ida||1999). In this method, the cutoff length corresponds to the box size L. 


2.2. Ring Model 


Ring particles are assumed to be identical for the sake of simplicity except for Section 4.4 
For the identical particle model, all particles have the same mass rrii = nip. This assumption is 


unrealistic since many observations show a power-law size distribution of ring particles (Zebker 


et al. 1985} French Nicholson 2000). The effect of the realistic size distribution is discussed in 
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Section 4.4 where we adopt the power-law size distribution (Zebker et al.|[l9^ French Sz Nicholson 


2000 ) 


dN 

dRn 


oc V’ 


( 6 ) 


where Rp is the particle radius and q is the power-law index. From the Voyager radio occultation 
measurements, the power-law index increases with the Saturnicentric distance from 2.7 to 3.0 in the 
A ring (Zebker et al.[1985). The stellar occultation observations showed the consistent results that 


the power-law indexes are q = 2.75 and 2.9 in the inner and outer A rings, respectively (French & 


Nicholson 2000). Considering these results, we choose the intermediate value q = 2.8. We set the 

= 150/^ cm and = 150^ cm where 7 = i?max/-Rmin = 


lower and upper size limits as R 
1 - 10 . 


For a ring with identical particles, the dynamical optical depth is given by 

__,,, / 7 ir\ \ / /T.v \ 

r = = 0.417 


^0 


Pp 


( 7 ) 


dRpPp \50gcm~2y ylOOcmy \0.9gcm~3 

where N is the number of particles in the simulation box and pp is the particle density. For a 


low optical depth ring, the dynamical optical depth is equal to the observed optical depth (Salo 


Sz Karjalainen 2003; Robbins et al. 2010), while for a high optical depth ring, it is larger than the 
observed optical depth. Though the particle density of solid ice is 0.9gcm“^, pp is estimated to 
be 0.45 g cm“^ from the best fit to the observations of azimuthal brightness variation in the A ring 
( jSalo et aT]|2004 French et akpOOT Porco et al.[2008 ). The surface density is inferred from density 
waves, which is typically Sq — 45gcm“^ in the A ring (e.g. Tiscareno et al. 2007). The surface 


density, radius, and particle density of ring particles collectively affect the structure of self-gravity 
wakes. We vary the dynamical optical depth from 0.3 to 1.8 in the identical particle models, which 
corresponds to the surface density of Sq = 36-216gcm“^ with Rp = 100cm and pp = 0.9gcm“^. 

Both the semi-major axis and particle density affect the spatial structure. If the distance of 
the simulation box from the planet, i.e., the semi-major axis is large, the tidal force becomes less 
effective, allowing ring particles to gravitationally accumulate. In the region outside the Roche limit, 
particles can form gravitationally-bound aggregates. The particle density affects the strength of 
the self-gravity. Larger density makes particles pull each other more strongly, causing more distinct 
wakes. These factors are combined into one non-dimensional parameter, the Hill radius of a particle 
pair relative to the sum of their physical radii, defined as 

-1/3 / „ \ 1/3 


ru = 


Rh 


= 0.82 


Ms 

5.69 X 1029 g 


Pp 

0.9 g cm“3 


a 

10 ^ km 


( 8 ) 


where Rh is the Hill radius of the particle Rh = (2mp/3Ms)^/^a. With the particle density pp 
fixed, th is proportional to a. We vary the ratio from 0.68 to 1.06. 


We introduce a parameter, initial Toomre Q to determine the initial random velocity of par¬ 
ticles, 

Qini — 


CT 


3.36GSn’ 


(9) 
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where Gx is the initial radial velocity dispersion (Toomre 1964). Since the ring relaxes into the 
stationary state quickly, Qini does not affect the final state (Salo 2012). Thus, we adopt Qini = 2. 

The ring models are summarized in Table for identical particles and in Table for particles 
with size distributions. Note that the dynamical properties of rings with the identical particles are 
characterized by the three non-dimensional quantities: r, fn and e (e.g., Karjalainen fc: Salo||2004 ). 
The dimensional quantities such as particle size are irrelevant to simulation results. However, in 
order to compare with the actual rings, we add the corresponding dimensional quantities to Table 
[^assuming = 100cm and pp = 0.45-0.90gcm“^. The surface density Sq is calculated by 
Equation Q and the Saturnicentric distance a is calculated by Equation Q. 









Table 1. Ring Parameters and Pitch Angles for Identical Particle Models 


Model 

T 

f-H 

e 

N 

So 

( gcm - 2 ) 

a 

( 10 ® km ) 

L 

( km ) 

di 

( deg .) 

( deg .) 

Xb 

(Acr) 

\ yh \ 

(Acr) 

1 

0.60 

0.81 

0.50 

44156 

36-72 

0 . 78 - 0.99 

0.48 

25.8 ± 4.1 

13.3 ± 4.2 

0.26 

0.53 

2 

0.70 

0.81 

0.10 

70119 

42-84 

0 . 78 - 0.99 

0.56 

27.3 ± 6.0 

10.9 ± 6.2 

0.28 

0.54 

3 

0.70 

0.81 

0.20 

70119 

42-84 

0 . 78 - 0.99 

0.56 

25.7 ± 4.3 

15.4 ± 7.0 

0.27 

0.55 

4 

0.70 

0.81 

0.30 

70119 

42-84 

0 . 78 - 0.99 

0.56 

26.0 ± 4.3 

13.8 ± 5.4 

0.31 

0.63 

5 

0.70 

0.81 

0.40 

70119 

42-84 

0 . 78 - 0.99 

0.56 

25.7 ± 3.9 

15.6 ± 4.7 

0.23 

0.49 

6 

0.70 

0.81 

0.50 

70119 

42-84 

0 . 78 - 0.99 

0.56 

25.2 ± 3.1 

15.1 ± 5.3 

0.26 

0.55 

7 

0.70 

0.81 

0.60 

70119 

42-84 

0 . 78 - 0.99 

0.56 

27.0 ± 3.1 

14.1 ± 3.3 

0.14 

0.28 

8 

0.70 

0.81 

< i ( Vn ) 

70119 

42-84 

0 . 78 - 0.99 

0.56 

24.5 ± 4.2 

13.9 ± 4.1 

0.27 

0.60 

9 

0.80 

0.60 

0.50 

17290 

48-96 

0 . 58 - 0.73 

0.26 

17.9 ± 1.5 

10.0 ± 1.5 

0.74 

2.29 

10 

1.00 

0.60 

0.50 

33770 

60-120 

0 . 58 - 0.73 

0.33 

22.1 ± 2.3 

13.5 ± 1.5 

0.04 

0.10 

11 

1.20 

0.60 

0.50 

58355 

72-144 

0 . 58 - 0.73 

0.39 

9 . 7 ± 2.1 

17.5 ± 6.1 

0.24 

1.39 

12 

1.40 

0.60 

0.50 

92666 

84-168 

0 . 58 - 0.73 

0.46 

8.2 ± 3.7 

19.8 ± 10.5 

0.19 

1.28 

13 

1.60 

0.60 

0.50 

138324 

96-192 

0 . 58 - 0.73 

0.52 

0.5 ± 2.7 

- 0.3 ± 2.3 

0.03 

3.98 

14 

1.80 

0.60 

0.50 

196950 

108-216 

0 . 58 - 0.73 

0.59 

9.8 ± 1.4 

17.8 ± 6.5 

0.29 

1.70 

15 

0.30 

0.81 

0.50 

5519 

18-36 

0 . 78 - 0.99 

0.24 

30.2 ± 3.0 

13.5 ± 2.6 

0.29 

0.50 

16 

0.40 

0.81 

0.50 

13083 

24-48 

0 . 78 - 0.99 

0.32 

30.3 ± 2.5 

13.2 ± 3.4 

0.22 

0.37 

17 

0.50 

0.81 

0.50 

25553 

30-60 

0 . 78 - 0.99 

0.40 

29.5 ± 2.8 

12.3 ± 4.4 

0.21 

0.37 

18 

0.70 

0.81 

0.50 

70119 

42-84 

0 . 78 - 0.99 

0.56 

25.2 ± 3.1 

15.1 ± 5.3 

0.26 

0.55 

19 

0.80 

0.81 

0.50 

104667 

48-96 

0 . 78 - 0.99 

0.64 

25.9 ± 3.9 

16.9 ± 6.0 

0.16 

0.32 

20 

0.90 

0.81 

0.50 

149028 

54-108 

0 . 78 - 0.99 

0.72 

24.5 ± 3.9 

12.8 ± 5.8 

0.35 

0.76 

21 

1.00 

0.81 

0.50 

204428 

60-120 

0 . 78 - 0.99 

0.80 

24.5 ± 4.8 

13.4 ± 5.7 

0.30 

0.67 

22 

1.10 

0.81 

0.50 

272094 

66-132 

0 . 78 - 0.99 

0.88 

23.4 ± 4.1 

13.6 ± 6.3 

0.47 

1.09 

23 

1.20 

0.81 

0.50 

353253 

72-144 

0 . 78 - 0.99 

0.96 

23.7 ± 3.9 

11.2 ± 6.2 

0.54 

1.23 

24 

1.30 

0.81 

0.50 

449130 

78-156 

0 . 78 - 0.99 

1.04 

24.7 ± 4.8 

13.6 ± 7.3 

0.37 

0.81 

25 

0.30 

0.68 

0.50 

1932 

18-36 

0 . 66 - 0.83 

0.14 

26.9 ± 5.2 

8.4 ± 5.2 

0.80 

1.58 

26 

0.30 

0.73 

0.50 

2957 

18-36 

0 . 71 - 0.89 

0.18 

26.4 ± 3.6 

13.5 ± 2.6 

0.44 

0.89 

27 

0.30 

0.77 

0.50 

4073 

18-36 

0 . 75 - 0.94 

0.21 

28.1 ± 2.9 

14.0 ± 3.3 

0.30 

0.57 

28 

0.30 

0.81 

0.50 

5519 

18-36 

0 . 78 - 0.99 

0.24 

30.2 ± 3.0 

13.5 ± 2.6 

0.29 

0.50 

29 

0.30 

0.85 

0.50 

7370 

18-36 

0 . 82 - 1.04 

0.28 

32.2 ± 2.7 

12.5 ± 3.2 

0.30 

0.47 

30 

0.30 

0.90 

0.50 

10386 

18-36 

0 . 87 - 1.10 

0.33 

37.1 ± 3.8 

13.1 ± 4.0 

0.27 

0.36 

31 

0.30 

0.94 

0.50 

13482 

18-36 

0 . 91 - 1.15 

0.38 

42.3 ± 4.0 

13.8 ± 4.1 

0.24 

0.26 

32 

0.30 

0.98 

0.50 

17312 

18-36 

0 . 95 - 1.20 

0.43 

47.3 ± 5.3 

13.3 ± 5.0 

0.24 

0.22 

33 

0.30 

1.02 

0.50 

22008 

18-36 

0 . 99 - 1.24 

0.48 

50.4 ± 6.0 

10.4 ± 5.2 

0.31 

0.25 

34 

0.30 

1.06 

0.50 

27722 

18-36 

1 . 03 - 1.29 

0.54 

53.7 ± 5.5 

7.3 ± 4.8 

0.41 

0.30 

35 

0.60 

0.68 

0.50 

15457 

36-72 

0 . 66 - 0.83 

0.28 

22.1 ± 1.3 

12.5 ± 1.6 

0.20 

0.50 

36 

0.60 

0.73 

0.50 

23660 

36-72 

0 . 71 - 0.89 

0.35 

26.4 ± 1.9 

13.3 ± 2.4 

0.11 

0.22 

37 

0.60 

0.77 

0.50 

32585 

36-72 

0 . 75 - 0.94 

0.41 

25.2 ± 2.3 

13.4 ± 4.0 

0.18 

0.39 

38 

0.60 

0.85 

0.50 

58965 

36-72 

0 . 82 - 1.04 

0.56 

27.9 ± 4.3 

13.7 ± 6.4 

0.32 

0.61 

39 

0.60 

0.90 

0.50 

83088 

36-72 

0 . 87 - 1.10 

0.66 

32.0 ± 5.6 

9.8 ± 6.3 

0.40 

0.64 

40 

0.60 

0.94 

0.50 

107858 

36-72 

0 . 91 - 1.15 

0.75 

35.7 ± 5.2 

16.7 ± 6.0 

0.25 

0.35 



3. 


Structure of Self-Gravity Wakes 


3.1. Formation of Self-Gravity Wakes 


Figure [^presents an example of the self-gravity wake formation for the standard model (model 
1). The parameters are r = 0.6, fn = 0.81, and e = 0.5. Initially particles are distributed uniformly 
in the simulation box (Figj^). 

Figure shows the time evolution of Q. The Q value is calculated from Equation ([^, where 
we calculate ax from the root mean square of the radial velocity dx/dt. Due to inelastic collisions 
between particles the kinetic energy decreases and the velocity dispersion and Q decrease. When Q 
becomes Q ~ 1-2, the non-axisymmetric perturbations grow significantly (Julian &: Toomre 1966 


Toomre 

1981 

). Accordingly, the condition for the emergence of the self-gravity wakes is obtained 

from Q <2 ( 

Salo 

1995 

Ohtsuki & Emori 

2000 

Daisaka et al. 

2001 

) 


> 


O.OSfjj^ for fn < 0.5 
0.2fH^/^ for rn > 0.5 


( 10 ) 


When this condition is satisfied, the self-gravity wakes appear finally regardless of Qini ( [S5^[m2| ). 
As shown in Fig{^, self-gravity wakes appear due to gravitational instability at i = 1.0 where 


t = t/(27r/D). Then the radial spacing of the self-gravity wakes is about Acr/2 (Salo et al. 2004) 


After reaching about unity Q increases with time and finally the time-averaged Q becomes 
almost constant ~ 2 after 5 orbital periods. In this state, the radial spacing is about Acr (Daisaka 


Sz Ida 1999 Salo et al. 2004). Since heating by gravitational scattering and cooling by inelastic 
collisions are balanced, the time-averaged velocity dispersion is constant. Figure [TJ: shows the 
equilibrium state of the self-gravity wakes. The high density region is trailing and inclined by 
10-30° with respect to the tangential direction. 

We performed simulations with various ring parameters and confirmed that the equilibrium 
state is reached after 5 orbital periods in all models if the axisymmetric overstability does not exist. 
Thus, hereafter we analyze the simulation result only for t = 5-20. 


3.2. Two-Dimensional Autocorrelation Function 


Since the self-gravity wakes are a transient and recurrent feature we need to investigate them 
in a statistical way. Two-dimensional autocorrelation function has been widely used to analyze 


self-gravity wakes (Toomre & Kalnajs 1991; Salo 1995; Daisaka & Ida 1999 Salo et al. 2004). We 


use the two-dimensional autocorrelation function of the surface density given by 

= / i:{x + x\y + y',t)i:{x\y',t)dx'(iy'-1. 

J-L/ 2 J-L /2 


( 11 ) 















































Table 1—Continued 


Model 

T 

fn 

e 

N 

So 

(gcm-2) 

a 

( 10 ® km) 

L 

(km) 

di 

(deg.) 

do 

(deg.) 

Xb 

(Acr) 

hbl 

(Acr) 

41 

0.60 

0.98 

0.50 

138497 

36-72 

0 . 95 - 1.20 

0.85 

39.0 ± 7.6 

14.0 ± 7.4 

0.31 

0.38 

42 

0.60 

1.02 

0.50 

176070 

36-72 

0 . 99 - 1.24 

0.96 

41.9 ± 7.2 

14.7 ± 6.8 

0.34 

0.37 

43 

0.60 

1.06 

0.50 

221779 

36-72 

1 . 03 - 1.29 

1.08 

46.4 ± 8.9 

6.7 ± 6.1 

0.56 

0.53 

44 

0.90 

0.68 

0.50 

52169 

54-108 

0 . 66 - 0.83 

0.43 

15 . 0 ± 1.7 

13.4 ± 1.8 

0.21 

0.78 

45 

0.90 

0.73 

0.50 

79854 

54-108 

0 . 71 - 0.89 

0.53 

21 . 6 ± 2.1 

14.8 ± 3.1 

0.09 

0.22 

46 

0.90 

0.77 

0.50 

109977 

54-108 

0 . 75 - 0.94 

0.62 

22.6 ± 3.2 

14.2 ± 5.3 

0.19 

0.45 

47 

0.90 

0.81 

0.50 

149028 

54-108 

0 . 78 - 0.99 

0.72 

24.5 ± 3.9 

12.8 ± 5.8 

0.35 

0.76 

48 

0.90 

0.85 

0.50 

199009 

54-108 

0 . 82 - 1.04 

0.83 

26.2 ± 6.1 

17.7 ± 7.1 

0.28 

0.57 

49 

0.90 

0.90 

0.50 

280423 

54-108 

0 . 87 - 1.10 

0.99 

31.2 ± 5.2 

14.7 ± 5.8 

0.30 

0.49 

50 

0.90 

0.94 

0.50 

364021 

54-108 

0 . 91 - 1.15 

1.13 

33.0 ± 9.9 

9.8 ± 7.6 

0.49 

0.76 

51 

0.90 

0.98 

0.50 

467429 

54-108 

0 . 95 - 1.20 

1.28 

39.0 ± 8.9 

14.2 ± 5.7 

0.28 

0.35 

52 

0.90 

1.02 

0.50 

594238 

54-108 

0 . 99 - 1.24 

1.44 

38.8 ± 11.6 

11.7 ± 9.4 

0.56 

0.70 

53 

0.90 

1.06 

0.50 

748505 

54-108 

1 . 03 - 1.29 

1.62 

47.4 ± 10.4 

14.3 ± 5.6 

0.29 

0.27 


(a) i = 0.0 (b) i = 1.0 (c) i = 5.0 
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Fig. 1.— Snapshots of the surface density normalized by the mean surface density Sq at t = 0.0 
(left), i = 1.0 (middle) and i = 5.0 (right) for model 1. 























Table 2. 


Ring Parameters and Pitch Angles for Size Distribution Models 


Model 

a 

(10® km) 

Pp 

(gem-®) 

-^min -Rmax 
(cm) 

N 

L 

(km) 

di 

(deg.) 

do 

(deg.) 

Xb 

(Acr) 

|Pbl 

(Acr) 

54 

1.30 

0.45 

150-150 

46148 

0.76 

30.0 ±5.0 

14.8 ±8.3 

0.21 

0.36 

55 

1.30 

0.45 

106-212 

47421 

0.76 

31.6 ±8.5 

15.6 ±6.9 

0.18 

0.30 

56 

1.30 

0.45 

87-260 

49392 

0.76 

34.3 ±8.0 

14.8 ±6.9 

0.20 

0.30 

57 

1.30 

0.45 

75-300 

52260 

0.76 

34.1 ± 8.0 

13.7 ± 6.4 

0.26 

0.39 

58 

1.30 

0.45 

67-335 

54745 

0.76 

35.5 ± 5.6 

13.4 ±5.3 

0.27 

0.38 

59 

1.30 

0.45 

61-367 

57128 

0.76 

37.5 ± 5.6 

14.9 ± 7.0 

0.23 

0.30 

60 

1.30 

0.45 

57-397 

58069 

0.76 

41.1 ± 8.5 

15.8 ±5.9 

0.21 

0.24 

61 

1.30 

0.45 

53-424 

60591 

0.76 

40.9 ±7.2 

14.7 ± 5.8 

0.24 

0.28 

62 

1.30 

0.45 

50-450 

62207 

0.76 

41.1 ± 6.9 

15.1 ± 6.1 

0.26 

0.29 

63 

1.30 

0.45 

47-474 

64929 

0.76 

41.6 ± 5.0 

15.1 ±4.7 

0.25 

0.28 

64 

0.85 

0.45 

47-474 

4998 

0.21 

22.0 ± 8.9 

18.2 ±5.5 

0.67 

1.66 

65 

0.90 

0.45 

47-474 

7043 

0.25 

28.0 ± 5.1 

11.0 ±4.2 

1.02 

1.92 

66 

0.95 

0.45 

47-474 

9742 

0.30 

22.6 ± 6.4 

18.8 ±3.9 

0.36 

0.87 

67 

1.00 

0.45 

47-474 

13254 

0.35 

34.2 ±4.7 

15.2 ±4.0 

0.26 

0.39 

68 

1.05 

0.45 

47-474 

17761 

0.40 

30.4 ± 5.7 

12.8 ±5.2 

0.41 

0.70 

69 

1.10 

0.45 

47-474 

23480 

0.46 

31.4 ± 4.3 

10.7 ± 5.7 

0.50 

0.82 

70 

1.15 

0.45 

47-474 

30657 

0.53 

32.8 ± 6.5 

13.7 ± 7.1 

0.30 

0.47 

71 

1.20 

0.45 

47-474 

39576 

0.60 

34.8 ± 5.3 

12.5 ±6.6 

0.32 

0.46 

72 

1.25 

0.45 

47-474 

50560 

0.68 

35.9 ± 7.1 

14.4 ±5.4 

0.29 

0.40 

73 

1.30 

0.45 

47-474 

63974 

0.76 

43.0 ± 5.7 

14.6 ±5.2 

0.25 

0.26 

74 

1.35 

0.45 

47-474 

80232 

0.85 

41.1 ± 8.2 

8.8 ±6.4 

0.48 

0.55 

75 

0.85 

0.90 

47-474 

2499 

0.21 

24.8 ± 13.9 

19.4 ±8.5 

0.89 

1.92 

76 

0.90 

0.90 

47-474 

3521 

0.25 

36.6 ± 11.0 

16.0 ±8.5 

0.53 

0.72 

77 

0.95 

0.90 

47-474 

4871 

0.30 

37.9 ±21.0 

14.8 ± 10.4 

0.34 

0.44 

78 

1.00 

0.90 

47-474 

6627 

0.35 

44.3 ± 16.2 

13.9 ±8.0 

0.48 

0.49 

79 

1.05 

0.90 

47-474 

8880 

0.40 

44.9 ± 13.4 

8.6 ± 7.9 

0.93 

0.93 

80 

1.10 

0.90 

47-474 

11740 

0.46 

44.3 ± 10.1 

15.4 ±9.2 

0.42 

0.43 

81 

1.15 

0.90 

47-474 

15328 

0.53 

45.6 ±26.1 

12.5 ±9.1 

0.49 

0.48 

82 

1.20 

0.90 

47-474 

19788 

0.60 

54.9 ± 7.0 

12.6 ±8.2 

0.48 

0.33 

83 

1.25 

0.90 

47-474 

25280 

0.68 

58.0 ± 8.0 

12.8 ±6.7 

0.51 

0.32 

84 

1.30 

0.90 

47-474 

31987 

0.76 

90.4 ± 8.8 

16.3 ± 11.2 

0.22 

0.00 

85 

1.35 

0.90 

47-474 

40116 

0.85 

90.2 ± 2.9 

16.2 ±8.7 

0.28 

0.00 


Note. — The surface density is So = 50 g cm ^ and the restitution coefficient is e = 0.5. 
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We calculate the autocorrelation with the uniform grid of 256 x 256. We first calculate the surface 
density at each grid from particle distribution and then evaluate the integral. Examples are shown in 
FigureThough there are particle aggregates with various shapes in the surface density snapshot, 
the autocorrelation analysis enables us to extract a smooth wake structure. As shown in Figure [^, 
b, and c, although the wake structure is almost steady, it slightly fluctuates with time. Therefore 
we introduce the time-averaged autocorrelation: 

^{x,y)= ^ ( 12 ) 

t2 — tl 

where we set = 5 and f 2 = 20. We approximate the time integral by the sum of 150 snapshots in 
i = 5.0-20.0 with At = 0.1. Figure shows the time-averaged autocorrelation over i = 5.0-20.0. 
The small variation is completely smoothed out. The inclined structure at the center corresponds 
to the structure of the high density region in the self-gravity wakes. This inclination means that 
the structure is trailing and the pitch angle is about 20°. The regions with negative autocorrelation 
coefficient correspond to voids between self-gravity wakes. 

It should be noted that the grid of 256 x 256 is reasonable. One might worry that the number 
of grids larger than that of particles leads to large density error. This is always true for low density 
parts where particles are sparsely distributed even if the number of particles are larger than that 
of grids. However, it is not problematic here since we focus on high density parts automatically 
considering the autocorrelation function. In addition, time-averaging of the autocorrelation function 
reduces the density error. To demonstrate a suitable grid size Figure 4 demonstrates the validity 
of the grid size. The function shape is correctly caught by the grid of 256 x 256. The typical 
discretization error of the pitch angle is smaller than ~ 4%. 


3.3. Calculation of Pitch Angle 


We calculate the pitch angle from the time-averaged autocorrelation. First, we extract the 
ridge of the dense part of autocorrelation by considering the integral curve of the gradient field 
of the autocorrelation function. Figure shows examples of the ridge of autocorrelation, which 
conhrms the previous results that the slope of the ridge line varies with the distance from the 
center ( Salo et ah]|2004 ). The slope in the outer part of the autocorrelation function is steeper than 
that in the inner part. Clearly we cannot fit the ridge by one line. Therefore we fit the ridge with 
the following function: 


y = — cot 6oX + (cot Oo — cot 0i)xberf ( ^ ) , 

V 2xb ) 


(13) 


where 0i, 0o, and Xb are parameters. The slope of the tangent to the fitting curve around the center 
\x\ <C Xb is cot^i and that for |x| ;:g> Xb is cot6*o. Hereafter, we define the inner and outer pitch 
angles as 0i and 9o, where 9i and 9o correspond to the pitch angles in the inner and outer parts of 
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Fig. 3.— Snapshots of the two-dimensional autocorrelation function at i = 5.0 (a), 10.0 (b), and 
15.0 (c) for model 1. The panel (d) shows the time-averaged autocorrelation over t = 5.0-20.0. 
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the autocorrelation function, respectively. The value Xb means the typical length where the pitch 
angle changes from 9\ to Oq. An example of the fitting function is shown in Figure 


We introduce a line integral of autocorrelation function along the ridge defined by Equation 
(13) as 

S’(6»i,6'o,Xb) = J ^{x,y)dl. (14) 

We calculate the integral by Simpson’s rule and find Oi, 9o, and Xb that maximize S by the downhill 
simplex method ( Nelder fc Mead||1965 Press 2007). 

We obtain the representative pitch angles from the time averaged autocorrelation function. 
However, the actual pitch angles fluctuate with time around the representative values. Thus, in 
order to estimate the fluctuation degree, we calculate the median absolute deviation from the 
time series of the inner and outer pitch angles calculated from the autocorrelation functions for 
t = 5.0-20.0. Examples of the optimized parameters and the median absolute deviation are 9i = 
(25.8 ± 4.1)°, 9o = (13.3 ± 4.2)°, and Xb/Acr = 0.26 for model 1 and 9; = (41.9 ± 7.2)°, 6 »o = 
(14.7 ± 6 . 8 )°, and x^/Xcr = 0.34 for model 42. As shown in Figure]^ the ridge curves agree well 


with Equation (13). Model 42 shows the deviation of the extracted ridge from the fitting curve 


in the low autocorrelation region. This is because that the ridge extracted by the gradient field 
is sensitively affected by the noise. In the low autocorrelation region, the htting curve is more 
robustly determined than the extracted ridge. 


The inner and outer pitch angles 9i and 9o, and the breakpoint position (xb, Uh) are summarized 
in Tables and where ?/b is calculated from yb = —Xh/ tan^i. The breakpoint Xb is generally less 
than Acr and about (0.1-0.5)Acr. 


4. Parameter Dependence 

We first adopt the identical particle model and examine the dependence of the pitch angle on 
the restitution coefficient e, the optical depth r, and the Hill radius of a particle pair relative to 
the sum of their physical radii fn. Next, we investigate the effect of the size distribution of ring 
particles. 


4.1. Restitution Coefficient 


Eirst we adopt the constant restitution coefficient model and investigate its effect on the pitch 
angle. The restitution coefficients adopted here are e = 0.1,0.2,0.3, 0.4, 0.5 and 0.6 (models 2-7). 
The other parameters are r = 0.7 and rn = 0.81. The results are summarized in Figure In all 
models, the inner and outer pitch angles barely depend on e, and their mean values are about 26.1° 
and 14.1°, respectively. As shown by Salo et al. (2004), the inner pitch angle is larger than the 


outer one. Since the influence of the self-gravity relative to the shear is strong around the center of 
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(a) 64 X 64 


(b) 128 X 128 


(c) 256 X 256 



Fig. 4.— Time-averaged autocorrelation functions for model 30 with grid of (a) 64 x 64, (b) 
128 X 128, and (c) 256 x 256. 



Fig. 5.— Autocorrelation function and the extracted ridge (dotted) for models 1 (left) and 42 
(right). The dashed curve shows the fitting function given by Equation (|13[) . 
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self-gravity wakes, the shape in the inner part of the autocorrelation function is nearly spherical. 
Conversely, the shape in the outer part is elongated due to the shear. Thus, we have 6\ > Bo¬ 
lt is expected that the pitch angle does not change for the velocity-dependent restitution 
coefficient model, since the pitch angles barely depend on e in the constant restitution coefficient 
models. This is confirmed by the simulation with the velocity-dependent restitution coefficient 
(model 8) where we assume a = 1.25 x 10^ km, pp = 0.45 gcm“^, and Rp = 120cm (Figj^. The 
inner and outer pitch angles are 24.5° and 13.9°, respectively. They are approximately the same 
as those for the constant restitution coefficient models. It is clearly shown that the pitch angle is 
independent of the restitution coefficient. 


4.2. Optical Depth 


We investigate the dependence of the pitch angle on the dynamical optical depth r with fixed 
th and e. The restitution coefficient is e = 0.5 and the Hill radius of a particle pair relative to the 
sum of their physical radii is fn = 0.6 and 0.81. The optical depths are r = 0.8-1.8 for fn = 0.6 
(models 9-14) and r = 0.3-1.3 for ru = 0.81 (models 15-24). The results are shown in Figure]^ For 
th = 0.6 the inner pitch angle roughly decreases with r, and becomes lower than 10°. For r > 1.0 
the outer pitch angle is larger than the inner one. However, it is arguable whether these values can 
be considered as the representative pitch angle of self-gravity wakes. Figure [^illustrates the surface 
density and its autocorrelation for r = 1.6 and fjj = 0.6 (model 13) at t = 19. We find that the 
axisymmetric structure and the non-axisymmetric self-gravity wakes coexist. This axisymmetric 


structure is caused by the overstable oscillation ( 

Schmit &; Tscharnuter 

[TM5 

Schmidt et al.||2001 

Salo et al.|2001 

Schmidt &: Salo|2003 

Latter & Ogilvie|2006 

). The axisymmetric structure appears 


for T > 1 and rn ~ 0.6. The overstable oscillation causes the small pitch angles ([French et al.j2007). 


On the other hand, for th = 0.81 the outer pitch angle does not depend on r, and for r > 0.5 
the inner pitch angle does not have clear dependence on r. For r > 0.5 the mean inner and 
outer pitch angles are 24.7° and 13.7°, respectively. For 0.3 < r < 0.5 the inner pitch angle is 
systematically larger than 24.7°. The inner pitch angles are 30.2°, 30.3°, and 29.5° for r = 0.3 
(model 15), 0.4 (model 16), and 0.5 (model 17), respectively. Figure 10 represents the surface 
density snapshots at t = 20 for r = 0.30 (model 15), 0.40 (model 16), 0.50 (model 17), and 0.60 
(model 18). From Equation (10), the condition for self-gravity wake formation is r > 0.27. Thus 
the ring of model 15 is marginally unstable and the self-gravity wakes are small and faint. As r 
increases (models 16, 17, and 18), the clearer long self-gravity wakes and voids develop. A possible 
explanation for the large inner pitch angle at small r is as follows. In marginally unstable rings, the 
self-gravity is weak and the resulting self-gravity wakes are faint. Since the wake clumps are loosely 
bound, they are easily destroyed by the interactions among clumps before they are elongated. In 
short we find that as long as the overstable oscillation does not occur and the clear self-gravity 
wakes form with large r, the inner pitch angle is almost independent of r. 
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40 
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30 
25 

th 

15 
10 
5 
0 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

6 

Fig. 7.— Inner and outer pitch angles 9i and Oq are plotted against the restitution coefficient e with 
r = 0.7 and rn = 0.81 (models 2-7). The open and filled symbols denote 9i and Oq, respectively. 
The solid and dashed lines denote 6i and 9o for the velocity-dependent restitution coefficient model, 
respectively (model 8). 
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Fig. 8.— Inner and outer pitch angles 9i (open) and 9o (filled) are plotted against the optical 
depth r. The restitution coefficient is e = 0.5 and Hill radius of a particle pair relative to the sum 
of their physical radii is fn = 0.6 (left, models 9-14) and 0.81 (right, models 15-24). 



Fig. 9.— Snapshots of the surface density (left) and its autocorrelation (right) for the model with 
T = 1.6, th = 0.6, and e = 0.5 (model 15) at i = 19.0. 
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X / Xgf X / 

Fig. 10.— Snapshots of the surface density normalized by the mean surface density at t = 20.0 
for models 15 (a), 16 (b), 17 (c), and 18 (d). 


(a) r = 0.30 (model 15) 


(b) r = 0.40 (model 16) 


XIX- 


(d) r = 0.60 (model 18) 


X ! 

(c) T = 0.50 (model 17) 
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4.3. Hill Radius Relative to Sum of Particles Radii 

The dependence of the pitch angle on fn is shown for r = 0.3 (models 25-34), 0.6 (models 


35-43), and 0.9 (models 44-53) in Figure 11 The outer pitch angle is independent of both r and 
fn. The mean outer pitch angles for r = 0.3, 0.6, and 0.9 are 12.0°, 12.8°, and 13.8°, respectively. 
This reason will be discussed in Section [5Tl 

On the other hand, the inner pitch angle increases with fjj from 30° to 50° for fn > 0.8. 
The similar tendency was shown in the previous numerical studies ( Salo|19^ Takeda fc Ida|[200T ). 
However, for fjj < 0.8 and r = 0.6 (models 35, 36, and 37), the inner pitch angle is almost constant, 
and for fjj = 0.77 and r = 0.3 (model 27) the inner pitch angle is unexpectedly large. The rings of 


these models are marginally unstable. Thus, as discussed in Section 4.2, the pitch angle becomes 
large. 


For 0.8 < fn 
pitch angle: 


using the least-square ht method, we obtain the fitting formula of the inner 


~ < 


(117.9 ± 12.0)rH 


(-65.1 ± 11.3) 

deg. 

(r = 0.3), 

(94.8 ± 4.9)fH 


(-52.0 ±4.6) 

deg. 

(r = 0.6), 

(113.3 ± 21.6)rH 


(-69.5 ±20.4) 

deg. 

(r = 0.9), 


(15) 


which shows that the dependence of 9\ on rn is independent of r. Thus, fitting the results of all r 
models we obtain a general formula for rjj > 0.8 and r = 0.3-0.9: 


0; ~ (108.7 ± 13.4)fH + (-62.2 ± 12.8) deg. 


(16) 


4.4. Size Distribution 


We investigate the effect of the size distribution of ring particles on the pitch angle. We adopt 
the power-law size distribution described by Equation Q with the lower and upper size limits 


t = 0.3 


T = 0.6 


t = 0.9 



Fig. 11.— Inner and outer pitch angles 9\ (open) and 9o (filled) plotted against rn with the optical 
depth r = 0.3 (left, models 25-34), r = 0.6 (middle, models 35-43), and r = 0.9 (right, models 
44-53). The solid lines denote the fitting formula given by Equations (16), and the dashed lines 
denote the estimation given by Equation (22) with = 0.45 and = 1.0. 
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(models 54-63). Figure 12 illustrates the dependence of the pitch angle on the ratio R m»^ /Rr, 


-2 


a = 


The surface density, the Saturnicentric distance, and the particle density are Sq = 50gem 
1.3 X lO^km, and pp = 0.45 gcm“^, respectively. The inner pitch angle increases with Rmax/Rmm- 
The inner pitch angle with the size distribution is 30.0°-41.6° and is 39% larger than that for 
identical particles. This is consistent with the previous results by Salo et al. (2004). As Rmax/Rmin 


increases, the filling factor and density of self-gravity wakes increase. The higher density of self¬ 
gravity wakes leads to the larger pitch angle because the shear in the higher density wakes is weaker 
relative to the self-gravity. The details of this dependence will be discussed in Section 5.2 On the 
other hand, the outer pitch angle barely depends on Rmax/Rmin- The mean outer pitch angle is 
14.8°. 


In Section 4.3 we found that the pitch angle for identical particles increases with fn- We expect 
the same trend for the pitch angle in the size distribution models. Figure [T^ shows the dependence 
of the pitch angle on the Saturnicentric distance. The surface density is Sq = 50gcm“^ and the 
particle density is pp = 0.45gcm“^ (models 63-73) and 0.9gcm“^ (models 74-84). The lower 
and upper size limits are fixes as i?min = 41.7 cm and i?max = 417 cm. The inner pitch angles for 
Pp = 0.45 gcm“^ are smaller than those for pp = 0.9gcm“^. This is because smaller pp corresponds 
to smaller fjj. As expected, the inner pitch angle generally increases with the Saturnicentric 
distance. For pp = 0.45 gcm“^ the pitch angle increases with the Saturnicentric distance a from 
22.0° to 43.0°, and the outer pitch angle has no trend and ranges within 8.8°-18.8°. Similarly, for 
Pp = 0.9gcm“^, the pitch angle increases with a from 24.8° to 58.0° for a < 1.25 x 10^ km, and the 
outer pitch angle is 8.6°-19.4°. 


14 


For Pp = 0.9gcm“° and a > 1.3 x 10° km, the inner pitch angle is about 90°. Figure 
presents the snapshot of the surface density and the autocorrelation for a = 1.35 x 10® km and 
Pp = 0.9gcm“^ (model 84). In the surface density snapshot, we cannot observe the distinct 
elongated wakes but nearly spherical aggregates. This is because the aggregate formation starts 
from a ~ 1.3 x 10® km for pp = 0.9gcm“^ 


Ohtsuki 

1993 

Salo||l995 

Karjalainen & Salo 

2004 


the autocorrelation function, the shape of the dense part is an ellipsoid, which corresponds to the 
Hill sphere. The long axis of the ellipsoid is almost parallel to the x-axis. Thus, in this case, the 
inner pitch angle becomes around 90°. The typical mass of the self-gravity wake is SqA^j.. Thus, 
supposing that a particle with mass SoA^ is located at the center, we obtain the corresponding 
Hill radius (i;oACT/3Ms)^/^a = Acr/(127r^)^/^ ~ 0.2Acr. This length scale is consistent with the size 
of the ellipsoid in autocorrelation function. 


We compare the simulation results with the observations. Hedman et al. (2007) found that 


the pitch angle increases from 17° to 26° with the Saturnicentric distance in the A ring. Colwell 


et al. (2007) found the similar results that the pitch angle in the A ring is 10-50° and increases 


with the Saturnicentric distance. Our simulations with pp = 0.45 gcm“® show that the inner pitch 
angle increases with the Saturnicentric distance around a ~ 1.3 x 10® km and the inner and outer 
pitch angles are 43.0° and 14.6°. We expect that in general the observed pitch angle ranges within 
the inner and outer pitch angles. Therefore, our simulations with pp = 0.45 gcm“^ are consistent 
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Fig. 12.— Inner and outer pitch angles 0^ (open) and 9o (filled) plotted against Rmax/Rmm 
for models 54-63. The surface density, the Saturnicentric distance and the particle density are 
So = 50gcm“^, a = 1.3 x lO^km and pp = 0.45 gcm“^, respectively. The solid and dashed lines 
denote the estimation given by Equation (22) with = 0.45 (solid) and = 1.0 (dashed), 
respectively. 
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Pp = 0.45 gem' 
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Fig. 13.— Inner and outer pitch angles 6i (open) and 9o (filled) plotted against the Saturnicentric 
distance. The particle densities are pp = 0.45 gcm“^ (left, models 63-73) and pp = 0.45 gcm“^ 
(right, models 74-84). The solid lines denote the estimation given by Equation (22) with Cw = 1.0 
and Cd = 1.0. 



Fig. 14.— Surface density snapshot at t = 20.0 (left) and its time-averaged autocorrelation (right) 
with the extracted ridge (dotted) and the fitting function (dashed) for model 84. 
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with the observations. 


5. Discussion 


5.1. Comparison with Galactic Spirals 


We compare the pitch angles of spiral arms in disk galaxies and self-gravity wakes in planetary 
rings. Disk galaxies are considered as collisionless systems. Recently, the pitch angle of spiral arms 


was investigated by the numerical and analytic studies (Grand et al. 2013, Michikoshi &: Kokubo 


2014). In the galactic disk, the pitch angle depends only on the shear rate defined by 


r = 2 - 


2D2^ 


(17) 


where k is the epicycle frequency. From the numerical simulations and the linear analysis [Michikoshi 


&: Kokubo (2014) obtained the fitting formula of the pitch angle 


7 r 


(18) 


For the shear rate of the Keplerian rotation F = 3/2, we obtain the pitch angle 9 ~ 11°, which is 
consistent with the outer pitch angle of self-gravity wakes. 


In the central dense parts where self-gravity wakes are closely packed, inelastic collisions play 
an important role in dynamics. Thus, the inner pitch angle of self-gravity wakes deviates from the 
estimated value for collisionless systems. In the outer low density regions, however, the effect of 
the shear is relatively strong compared to that of collisions, and the shear rate mainly determines 
the outer pitch angle just like collisionless systems. Thus, the outer pitch angle is the same as that 
of collisionless systems and is independent of any ring parameters. 


5.2. Inner Pitch Angle 

We simplify the physical processes concerning the inner pitch angle and derive a rough estimate 
as follows. In high density regions, the self-gravity is strong. Without the shear and only with the 
self-gravity, an ellipsoidal aggregate whose shape is the Hill sphere would form and its long axis is 
parallel to x-axis, that is, the pitch angle is 90°. Therefore, the strong self-gravity tends to increase 
the inner pitch angle to 90°. Conversely, the shear rotates the aggregate and tends to decrease the 
inner pitch angle to 0°. The inner pitch angle may be determined by the balance between these 
two processes. 

We consider the evolution of an aggregate formed due to the self-gravity. When the aggregate 
is deformed by the shear, the self-gravity tends to restore it back to its original shape. The 
restoration timescale may be approximated by tgrav ~ 1 / \/Gp„ where is the bulk density of 
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self-gravity wakes, which is the timescale of the Jeans instability or gravitational instability (e.g., 

. The bulk density pw is less than the particle density pp. Introducing 
a dimensionless factor C*, we assume pw = CwPp- For identical particles, C* is less than 0.74 that 
is for the closest packing. The timescale of self-gravity is given by 


Chandrasekhar &: Fermi| 1953) 


f 

^grav — 


Cd ^ Q 



TT _ 


3/2 


(19) 


where Cd is a non-dimensional factor on the order of unity. The shear rotates and extends the wake 
and the pitch angle decreases with t as 


tan0 = 


3QC 


( 20 ) 


where t is the time elapsed since 9 = 90° (e.g., Toomre||1981 Binney Sz Tremaine||2008 ). The shear 
timescale is 


^shear 


1 d tan 6 


tan 6 dt 


-1 


30 tan 0 


As the pitch angle decreases with time, the wake rotation slows down. 

When the wake pitch angle is large (9 ~ 90°), the wake rotation is fast (tshear < 


grav; 


( 21 ) 


The 


wake rotates quickly before the self-gravity restores the inclined wake back to its original shape. As 
the pitch angle decreases, the wake rotation slows down and the self-gravity becomes important. 
Thus, from tghear = igrav we obtain a crude estimate of the inner pitch angle: 


tan^i ~ = 0.84 

CdV vr H 




Cd \0.9gcm~3 


1/2 


a n 3/2 
10® km/ 


( 22 ) 


which means that the pitch angle increases with C^. The bulk density of self-gravity wakes with 
size distribution can be larger than that for identical particles, and thus the pitch angle with 
size distribution is larger than that for identical particles, which is consistent with the simulation 
results. Figure 12 presents the pitch angle as a function of Rm&xlRm\n- As Rmax/Rmm increases, 
the hlling factor increases. Assuming Cd = 1.0, we hnd that = 0.45 and 1 are suitable 
for the identical particle models and the size distribution models of 7?max/7?min = 10, respectively. 


As shown in Figures 11 and 13 Equation (22) can roughly explain the dependence of the inner 


pitch angle on fn and a with C* = 0.45 for the identical particle models and = 1.0 for the size 
distribution models with J?max/7?min = 10. The simulation results and Equation (22) show that the 
inner pitch angle increases with fjj. Since rn is proportional to the Saturnicentric distance a, as fu 
increases, the orbital period increases, that is, the shear timescale increases. On the other hand, the 
timescale of the gravitational instability depends only on and is independent of a. Therefore, 
as fn increases, the self-gravity dominates over the shear and the inner pitch angle increases. 


Note that the actual hlling factor has not been investigated in detail yet. These results do 
not suggest that the hlling factor is = 0.45 or 1.0. The factor is also a free parameter. If 
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we adopt smaller Cd, the corresponding filling factor becomes smaller. The realistic filling factor 
should be smaller than unity. 


We used Rmax/Rmm < 10 in the size distribution models. However, this size ratio is far 
smaller than the observationally inferred values that are 70 and 2000 in the inner and outer A 
rings, respectively (French &: Nicholson 2000). As shown in Figure 12, the inner pitch angle is 
almost constant for R m»^ lRm\n > 5. Accordingly, the simulation results with R m»^ /R m\u = 10 in 
Section 4.4 may be applicable to the realistic ring. 


6. Summary 


We have investigated the pitch angle of the self-gravity wakes in the dense planetary ring by 
performing local Wbody simulations where both inelastic collisions and self-gravitational inter¬ 
actions between particles are included. In order to extract the feature of the self-gravity wakes, 
we calculated the two-dimensional autocorrelation function of the surface density. We measured 
the angle of the inclined structure of the autocorrelation function as the wake pitch angle. We 
obtained two pitch angles in the inner and outer parts of the autocorrelation function separately. 
We conhrmed that the inner pitch angle is larger than the outer one in all models. 


We examined the dependence of the pitch angle on various ring parameters. We found that the 
outer pitch angle is independent of the ring parameters, which is 10°-15°. This value is consistent 
with the estimation of pitch angle of the spiral arms in collisionless systems (Julian &: Toomre|1966 


Michikoshi & Kokubo 2014). The inner pitch angle does not depend on the restitution coefficient 


e of particles and the optical depth r provided that the overstable oscillation does not occur. We 
found that the inner pitch angle is determined by the Hill radius of a particle pair relative to the sum 
of their physical radii fn- The inner pitch angle increases with th from 20° to 50°. The tendency 
that the pitch angle increases with fu is consistent with the result of the observational studies 
that the pitch angle increases with the Saturnicentric distance (Hedman et al.|[2007 Colwell et al. 


2007). The pitch angle with the size distribution model is larger than that for the identical particles 


because the filling factor for the size distribution model is larger and the resulting self-gravity is 
stronger. The strong self-gravity increases the inner pitch angle. 


Considering the effects of the self-gravity and the shear, the tendency of the increase of the 
inner pitch angle with the Saturnicentric distance can be explained. However, this model is far 
from a complete theory. We need a more sophisticated theoretical framework to understand the 
overall activity of the self-gravity wakes and more detailed analyses of the self-gravity wakes by N- 
body simulation. It is especially important to understand the wake amplification and destruction 
processes considering the inelastic collisions and the nonlinear effects. We will focus on these 
processes in the next paper. 


Numerical computations were carried out on the GRAPE system at Center for Computational 
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